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This publication will deal with an explicit determination of the time evolution of the spin orienta- 
tion axes and the evolution of the orbital phase in the case of circular orbits under next-to-leading 
order spin-orbit interactions. We modify the method of Schneider and Cui proposed in ["Theo- 
reme liber Bewegungsintegrale und ihre Anwendungen in Bahntheorien" , Verlag der Bayerischen 
Akademie der Wissenschaften, volume 212, 2005.] to iteratively remove oscillatory terms in the 
equations of motion for different masses that were not present in the case of equal masses. Our 
smallness parameter is chosen to be the difference of the symmetric mass ratio to the value 1/4. 
Before the first Lie transformation, the set of conserved quantities consists of the total angular 
momentum J, the amplitudes of the orbital angular momentum and of the spins, L, Si, and S2. In 
contrary, 5 1 := | Si + S2 1 is not conserved and we wish to shift its non-conservation to higher orders 
of the smallness parameter. 

We perform the iterations explicitly to first order, while performing higher orders would mean 
no structural difference or harder mathematical difficulties. To apply this method, we develop a 
canonical system of spin variables reduced by the conservation law of total angular momentum, 
which is imposed on the phase space as a constraint. The result is an asymptotic series in e that 
may be truncated appropriately considering the physical properties of the regarded system. 
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I. INTRODUCTION 

Compact binaries are celestial systems that are likely 
to possess spin. Astronomical observations even sug- 
gest that accreting black holes can be spun up to more 
than 98% of the maximal (extremal) spin [1]. Further, 
compact objects are important sources for gravitational 
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waves, and also important for relativistic astrophysics 
[2]. A successful analytical method to deal with com- 
pact binaries in general relativity is the post-Newtonian 
(PN) approximation, which is applicable when the dis- 
tances are large and the velocities small compared to the 
speed of light c. The PN framework reduces the com- 
plicated Einstein equations (nonlinear partial differential 
equations) to ordinary differential equations. These PN 
equations of motion are often encoded in the form of a 
Lagrangian potential or a Hamiltonian. The goal of the 
present paper is to find approximate solutions to these 
equations of motion in the case of circular orbits when 
the spins of the objects are involved. The method em- 
ployed makes crucial use of a phase space structure, so wc 
employ PN results derived within the ADM (=Arnowitt- 
Deser-Misner) canonical formalism of general relativity, 
whose direct outcome is a description of the spin motion 
and the motion of the point mass with positions x Q and 
momenta p in terms of a Hamiltonian. For the ADM 
approach to the PN approximation in the presence of spin 
see [3, 4] and references therein. 

First attempts on spin in the PN approximation have 
been done in [5-8], where, e.g., the leading-order (LO) re- 
sults for spin-orbit and spin(l)-spin(2) effects have been 
calculated. Most important for the present paper are 
spin-orbit interaction Hamiltonians up to and including 
next-to-leading order (NLO), which were derived within 
the ADM formalism in [9, 10]. Corresponding equations 
of motion were already known before [11, 12], and other 
approaches succeeded at this order, too [13-15]. A gen- 
eralization of the NLO spin-orbit Hamiltonian to arbi- 
trary many objects was calculated in [16]. Even the next- 
to- next-to-leading order (NNLO) spin-orbit Hamiltonian 
was derived [17] and was recently confirmed [18, 19]. Sim- 
ilar results also exist for spin(l)-spin(2) interaction. The 
complete NLO Hamiltonian was calculated in [20] and 
confirmed by [21, 22] (a partial result is contained in 
[23], see [20, 21] for a discussion). The NNLO spin(l)- 
spin(2) interaction was simultaneously derived in a re- 
duced Hamiltonian form [24] and in a potential form [25] 
(a comparison is still missing, see [26] for the emerging 
difficulties). For the LO spin(l)-spin(l) interaction, see, 
e.g., [27]. The extension to NLO succeeded in the form 
of a fully reduced Hamiltonian [28-30] and in the form of 
a potential [31, 32]. The LO Hamiltonians of cubic and 
quartic order in spin were derived in [29, 33] (though 
crosschecks in [34] suggest that the quartic order is not 
complete yet). 

The determination of the far-zone gravitational waves 
generated by compact binaries requires the knowledge 
of certain radiation multipoles. The LO spin-orbit and 
spin(l)-spin(2) contributions were derived in [35]. Based 
on the energy and angular momentum flux obtained 
there, dissipative effects on the orbital evolution due to 
spin were inferred [36] (see [37] for the spin(l)-spin(l) 
level). Corresponding PN equations of motion can also 
be calculated directly by solving the Einstein equations in 
the near-zone [38, 39], see also [40, 41] for dissipative (ex- 



plicitly time-dependent) Hamiltonians. The NLO spin- 
orbit effects on the energy flux and the consequences for 
the evolution of the phase were derived in [42]. Source 
multipole moments up to quadratic order in the spins 
sufficient for the 3PN energy flux [43] and the 2.5PN 
radiation field [44] were calculated. The spin contribu- 
tions to the gravitational wave form were extended to 
2PN order explicitly in [45]. Even some spin-dependent 
hereditary contributions at 3PN were derived for circular 
orbits [46]. 

There already exist extensive banks of wave forms for 
circular binaries without spin, for example to be found 
in [47]. Gravitational waves from eccentric compact bi- 
naries have been calculated in [48] (see also references 
therein) to 3PN accuracy in the orbital motion. Re- 
cently, analytical gravitational wave form expressions in 
the time-Fourier domain for non-spinning coalescing bi- 
naries have been calculated through 2PN in [49, 50]. As 
one includes spin, the treatment of the orbital motion 
gets more complicated, because in general, spin preces- 
sion equations for the orbital angular momentum L have 
to be taken into account. For the case that the spins 
are aligned (up-up, down-down, or up-down configura- 
tions), the binary motion and the GW expressions are 
given in [51]. As soon as the spins are not aligned to 
L, the spins start precessing in a complicated manner 
such that until now, PN-exact solutions for the spin mo- 
tion at leading-order spin-orbit interaction are known for 
eccentric orbits only for two cases: (1) the single-spin 
case, and (2) the equal-mass case [52]. Very recently, in 
[46] an approximative solution to the circular-orbit prob- 
lem through leading-order spin-orbit interaction has been 
solved, neglecting specific higher-order terms emanating 
in the time evolution of their precession angles as well. 
We will extend this through next-to-leading order spin- 
orbit interaction as they follow from recent developments 
in the ADM formalism. Thus, we work with canonical 
variables throughout this article. 

The equations of motion following from the men- 
tioned Hamiltonians are ordinary differential equations 
and solving them numerically is straightforward and fast. 
However, for the purpose of data analysis this may still 
be too slow if overlap integrals with theoretical wave- 
forms for many initial conditions need to be calculated, 
e.g., [53]. In particular the presence of spins drastically 
increases the dimension of the parameter space and thus 
the number of initial conditions that must be considered. 
Therefore analytic solutions to the equations of motion 
including spin [52] are of great importance. Approximate 
analytic solutions can also be more accurate than the use 
of standard numerical techniques, especially when one 
evolves the system over several thousands of orbits. Such 
situations are not only relevant for gravitational wave as- 
tronomy, but also for predicting the distribution of recoil 
( "kick" ) velocities for merging black hole binaries due to 
the emission of gravitational waves [54] . (There the evo- 
lution of the spin orientations during the long inspiral 
period plays a crucial role.) Analytic solutions are also 
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important for hereditary effects (e.g., tails), as these de- 
pend on the full evolution history. The solutions derived 
in the present paper can therefore be useful to extend the 
spin-dependent tail effects (for circular orbits) presented 
in [46] to higher orders. As it seems to be impossible 
to cover the whole parameter space analytically, further 
approximations are necessary, such as expansions around 
the equal-mass [55] or circular-orbit 1 cases. 

As long as no radiation reaction effects due to the emis- 
sion of gravitational waves are regarded, the total angu- 
lar momentum defined as the sum of the individual spins 
and the orbital angular momentum L is conserved. This 
can play the role of a "weak" constraint in the sense of 
Dirac [56], also see [57]. The goal is to find a new set 
of generalized coordinates that fulfill standard canonical 
Poisson brackets. This could be done with the help of 
Dirac brackets. But if one is able to find generalized 
canonical coordinates that can eliminate the constraints 
by construction, one is not forced to walk this way. The 
hope is that, by reducing the number of dynamical vari- 
ables with the help of conservation laws, the problem of 
finding solutions to the equations of motion simplifies. 

Note that the magnitude of L will not be conserved 
if spin(l)-spin(2) interactions are included. We conclude 
that, for those terms, the Hamiltonian must depend on 
the conjugate coordinate to L, and the full dynamics is 
not described by precessions only. Anyway, as we include 
purely spin-orbit terms, the problem simplifies drastically 
and collapses to the dynamics of a sequence of precessions 
for a first insight. 

The present paper is organized as follows. After in- 
troducing dimensionless quantities and some notation in 
section II, we present the details of the construction of 
a reduced phase space for two objects with spin in Sec- 
tion III. We find that the three spin angle coordinates 
used in Ref. [55] are canonical, and their corresponding 
conjugated momenta are given by the total spin ampli- 
tudes along their rotation axes. Next, we show details of 
the Lie transformation algorithm in order to shift peri- 
odic terms to higher orders of the smallness parameter 
(here: deviation from the equal-mass case) in a pertur- 
bation theory in Section IV. This method is a modified 
version of the "modified Cui method" [58, 59]. This pro- 
cedure is applied to the binary spin-orbit Hamiltonian at 
next-to-leading order in Section V. There we also apply 
the method of shifting perturbative contributions to the 
Hamiltonian to even higher orders of the smallness pa- 
rameter. A very short review how to combine multiple 
Lie transformations is provided in Section VI, for those 
readers who like to go further than us. Final conclusions 
and outlook are given in VIII. 



II. DIMENSIONLESS QUANTITIES 

For a binary system it is suitable to work with dimen- 
sionless quantities only, which will be used throughout 
the paper. This can be achieved by measuring masses in 
terms of the reduced mass p of the binary, distances in 
terms of GM/c 2 , and time intervals in terms of GM/c 3 . 
Here M is the total mass of the binary, G is the Newton 
gravitational constant, and c is the speed of light. All 
quantities appearing in the present paper are assumed to 
be measured in these units from now on and are there- 
fore dimensionless. The only exception is the dimension- 
less Kerr parameter of black holes, which is defined by 
Xa = cSa/Gm 2 where S a is the spin-length and m a the 
mass of the a-th object. This relation defines Xa also for 
objects other than black holes, e.g., neutron stars. As- 
trophysically relevant cases are covered by Xa S 1- We 
therefore talk of rapid rotation if % a ~ 1 and of slow 
rotation if Xa -C 1. 

The dimensionless masses of the binary's constituents 
may be written as 



mi = 1 + p , m 2 = 1 



(2.1) 



where p = ~ is the mass ratio. Obviously all expressions 
can be written in terms of the mass ratio p only instead 
of the individual masses mi and m 2 . We also make use 
of the symmetric mass ratio 



>1 : 



171x1712 



P 



(mi+m 2 ) 2 (l + pY 



(2.2) 



which takes on the value | in the equal-mass case and 
tends to zero if one of the masses is much smaller than 
the other. We further introduce a parameter e by 



2 1 



(2.3) 



which measures the deviation from the equal-mass case. 
It holds 



1 - 2e 
l + 2e 



mi 



1 + 2e 



, m 2 



1 - 2e 



(2.4) 



where we assumed mi < m 2 . 

The relation between dimensionless spin-lengths S a 
and dimensionless Kerr parameters Xa of black holes 
reads 



S\= p Xi , S 2 = PX2 ■ 



(2.5) 



Manuscript unpublished, under preparation. 



Notice that the dimensionless spins S a are defined anal- 
ogous to the Xa m some other publications, e.g., in [55]. 

For the sake of simplicity and to avoid introducing a 
new smallness parameter, c is regarded to be the book- 
keeping parameter for the post-Newtonian approxima- 
tion and does not have the dimension of speed in the 
Hamiltonians for the remainder of this article. 
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III. CANONICAL VARIABLES 

In this section we introduce the canonical variables, or 
the phase space, of the binary system. This phase space 
structure is prerequisite to the application of Lie-series 
techniques for solving the equations of motion. In order 
to simplify the problem we reduce the dimension of the 
phase space by making use of certain conservation laws. 
This reduction is most transparent in a specific basis on 
phase space, which is derived in the present section. 

A. Center of Mass System 

In this section we illustrate how the number of phase 
space variables can be reduced by making use of conser- 
vation laws. We discuss this reduction by looking at the 
well known transition to the center-of-mass frame. The 
spins are neglected for now. 



Poisson brackets can be "read off" from the kinematic 
term, i.e., 

{r i ,p j }=5 ij , (3.5) 

all other zero. 

The phase space was reduced from 12 to 6 variables 
with the help of 6 conserved quantities (total linear mo- 
mentum and center-of-mass position) . Formally the con- 
servation laws can be treated as constraints on the phase 
space. However, such constraints are already fulfilled by 
the equations of motion, i.e., they do not produce ad- 
ditional constraint forces. Instead of working with an 
action, constraints on a phase space can be handled us- 
ing the Dirac bracket [56, 57, 60]. In this sense (3.5) can 
be considered as a Dirac bracket. 



2. Orbital angular momentum 



1. Reduction 

Consider an action W of the form 

W = J dt(p 1 -z 1 + p 2 -z 2 -H), (3.1) 

with particle positions z a and canonical momenta p a . 
The equations of motion follow from a variation of the ac- 
tion and take on the form of Hamilton's equations, where 
H plays the role of the Hamiltonian. (Remember that the 
variations of z a and p a are considered as independent.) 
The Poisson brackets read 

{z\, Plj } = Sij , {z l 2 ,P2j} = Sij , (3.2) 

all other zero. We call the terms in the action involving 
time derivatives kinematic terms in the present paper, as 
they correspond to the Poisson bracket structure. The 
interaction is described by the Hamiltonian. 

The dimension of the phase space is 12. This dimen- 
sion can be reduced by using conserved quantities. First 
the total linear momentum P is constant (to the post- 
Newtonian order considered here). The center-of-mass 
system is defined by a vanishing total linear momentum, 

P = Pi+p 2 = 0. (3.3) 

Then the position of the center-of-mass is also constant 
and may be chosen to be the coordinate origin. The 
action now reads 

W = J dt (p • r - H) , (3.4) 

where p = Pi = — P2 and r = zi — z 2 . Due to transla- 
tion invariance the Hamiltonian H depends on Zi and z 2 
solely through the combination r = z x — z 2 . The phase 
space now consists of the six variables r and p only. The 



We can reformulate the phase space by introducing an- 
gles for r and (orbital) angular momenta. We introduce 
a rotation matrix Ajq, 

transforming between the frame O, which co-rotates with 
the orbital motion of the binary, and the inertial frame 
/. Notice that Ay o can be parametrized by three angle 
variables, e.g., the Euler angles. We can choose the 1- 
axis of the co- rotating frame A]q such that it points in 
the direction of r, 

Ko = n l , (3.7) 

where n = r/r and r = |r|. The kinematic terms then 
turn into 

p • r = (n • p)r + y l k%K%kf r , (3.8) 
= p r f+~L ki n%, (3.9) 

where 

p r :=n-p, L ki :=2r^, : = A ZO A ZO • ( 3 - 10 ) 

Notice that the angular velocity tensor fijQ is antisym- 
metric, Qjq = — ^/o- L^i is the orbital angular momen- 
tum tensor. Corresponding vectors are given by 

L l := ^ijkLjk, := 2 e ijk^io ■ (3- 11 ) 

The implications of this new form of the kinematic terms 
for the phase space structure are discussed in the follow- 
ing. 
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B. Canonical angles: general theory 

In the center-of-mass frame, most of the dynamics of 
a binary system with spinning components can be de- 
scribed in terms of angular momenta, namely the orbital 
angular momentum and the two spins of the components. 
Before we derive canonical variables for this system, let 
us prepare some general developments on the theory of 
classical angular momenta. 



Angular momentum algebra 



We consider a generic angular momentum represented 
by an antisymmetric tensor Sij = ~Sji and a rotation 
matrix Ay c transforming to a "co-rotating frame C. In- 
spired by the kinematic terms found in the last section, 
we consider an action of the form 



W 



dt 



2 l 3 ic 



H(Ay c ,Si 



(3.12) 



Let us derive the equations of motion by independently 
varying Sij and the angle variables parameterizing Ay c . 
The details of this process are analogous to the relativis- 
tic case discussed in [60] and are not repeated here. (For 
example, it is easiest to use the independent antisymmet- 
ric variation symbol 89 % i = —89 : > l = A^ C #A^ C instead of 
explicitly varying the angle variables.) Again, the equa- 
tions of motion can be written as Hamilton's equations 
(with Hamiltonian H ) if we impose the Poisson brackets 

{A%, A%} = , {A%, S kl ] = 8 a A% - 8 ik A% , 

(3.13) 

{S^, Ski} = 8ikSji — SjkSu — SuSjk + SjiSik ■ (3.14) 

Notice that the last relation is the usual angular momen- 
tum algebra, 

{S i ,S*} = e iJk S k , (3.15) 

where S l :— h e ijkSjk is the spin vector. If the Hamil- 
tonian H does not depend on the orientation of the co- 
rotating frame, or A l T 3 c , then only this angular momen- 
tum algebra is needed. However, for now we look at the 
most general case where the phase space consists of three 
pairs of canonical variables contained in A l y c and Sij. 



2. Euler angles 

Before we can make the phase space structure more 
explicit, we need a way to parametrize a generic rotation 
matrix. One possibility is 

A 13 (oj, n) = n l n 3 + (5« — ri l n J ) cos a; — €ijkn k s'muj , 



which gives a rotation of angle u> around a unit vector 
n, which intrinsically provides the rotation axis. How- 
ever, throughout this paper we will parametrize rotation 
matrices in terms of three Euler angles a, /?, and 7 by 

A«(a,/3, 7 ) = A lfc (a,e 3 )A H (/3,e 1 )A^( 7 ,e 3 ), (3.17) 

where e a = (<5*). Notice that it holds A oi (a, — 
A lJ ; (— 7, — /3, —a) . The Euler angles always enter the kine- 
matic terms through the angular velocity 



<>'(n. I. V) = ie Ijfc A^(a,/3,7)A H (a,/3, 7 ). (3.18) 



A straightforward calculation leads to 

(cos a $ + sin a sin (3 7 
sin a p - cos a sin £7 | • (3.19) 
a + cos j3 7 

This result will be used frequently throughout the next 
sections. 



3. Canonical angle variables 

Let us introduce another frame S in which the 3-axis is 
aligned to Si (which in general is not necessarily identical 
to the co-rotating frame, e.g., for a non-spherical top). 
We parametrize the transformation as 



such that 



sin 9s sin 0s 
(S r ) = S(Ai 3 s ) = S [ - sin 9 S cos <j> s 

cos 9 S 



(3.20) 



(3.21) 



where S := |S|. We can then decompose the transforma- 
tion to the co-rotating frame C as Ay c = A^gA^ and 
parametrize A^ in terms of Euler angles, 



A^ c = A^(a s ,/3 s , 7s ). 



(3.22) 



(3.16) 



However, as the first and the last of the Euler angles is a 
rotation around the 3-axis, it holds 

A tk ^s,9s,0)A k ^a s ,/3s, ls ) 

= A 4fc (0 s , 9 S , a s )A k i(0, p s , 7s ). (3.23) 

Thus an equivalent formulation is 

A% = A^s, 9 S , a s ) , Ay c = A«(0, S> 7s) , (3.24) 

which in hindsight of later calculations is used from now 
on. The reason is that if the rotation axis of the co- 
rotating frame is parallel to the spin, then fts and 7s are 
constant and can usually even be set to zero. Then with 
the convention (3.24) the frames S and C coincide. 
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Using the decomposition A S 3 C = A IS Aj° c we obtain an 
"addition theorem" for the angular velocity, 

^sc = ^W^-is^is + ^si ' (3.25) 

with the angular velocity of the frame C relative to 
the frame S, 0g C := Ag C A^ c , and the angular ve- 
locity of the inertial frame / relative to the frame S, 
■= A^A^J (where A^j = Ajg). In terms of angular 
velocity vectors this reads 

n sc = A%n{ c + n* sl , (3.26) 

where we used £jjfcAj S A" l s A" s = ei mn . We can now 
rewrite the kinematic terms in the frame S, 

-Sij£l IC — S Sl IC — SA IS Vl l IC , (3.27) 
= Sfl 3 sc - Sn si . (3.28) 

Using 

n SI = w(-a s , -e s , -4> s ) , n sc = n*(o, p s , is) , 

(3.29) 

and (3.19) we finally have 

\*.M'u = S cos sis + Sa s + S cos 6 s 4>s , (3.30) 

= S 3 <p s + Shis + Sa s , (3.31) 

where S 3 := S cos 9s is the 3-component of the spin vec- 
tor in the inertial frame and Sq := S cos /3s is the 3- 
componcnt of the spin vector in the co-rotating frame. 
The Poisson brackets read 

1 = {0 S , S 3 } = { ls , S 3 C } = {a s , S} , (3.32) 

all other zero. We have explicitly expressed Ay c = 

A^gAg^ and Sij in terms of three pairs of canonical vari- 
ables. 

In the case that the Hamiltonian H is independent 
of the orientation of the co-rotating frame, that is, of 
75, Sq, and as, the phase space can be further simpli- 
fied. The canonical conjugates Sq, 75, and S are then 
cyclic and therefore constant. The only relevant Poisson 
bracket is {<f>s, S 3 } — 1 in this case. 

4- Example: Orbital angular momentum 

As an example we apply our findings from the last sec- 
tion to the orbital angular momentum. We parametrize 
the frame L, in which the orbital angular momentum is 
parallel to the 3-direction, as 

A i I i=A ij {<f> L) e L ,a L ) ) (3.33) 

such that 

(sin 9 L sin 4> L \ 
- siii 0£ cos 0l . (3.34) 
cos L J 



As L is orthogonal to n, we can choose the angle «l such 
that n is aligned with the 1-axis in the L frame. Then 
the frame co-rotating with the orbit O can be chosen to 
be identical to the frame L, A l L J = A u (0,0,0), as this 
choice satisfies the condition (3.7). According to the last 
section it holds 

h H nf Q = L$l\ L = L 3 j> L + La L , (3.35) 

and the Poisson brackets read 

l = { ( f >L ,L s } = {a L ,L}, (3.36) 

all other zero. 



5. Addition of angular momenta 

We now consider the case that we have two spins Si 
and S 2 which add up to a total spin S := Si + S 2 . Of 
course, we could just construct the phase space by two 
copies of (3.32). However, we try to include the total 
spin S in the phase space here. 

We utilize frames denoted by Si, S2, and S where 
the corresponding angular momenta point into the 3- 
direction, respectively. As the vectors Si, S2, and — S 
form a triangle and thus lie in a plane, it makes sense to 
introduce a frame in which this plane is fixed as, say, the 
2-3-plane. This can be achieved by rotating the frame S 
around S by a suitable angle as, so we again have 

A%=A*^<f> 8 ,e 8 ,a 8 ). (3.37) 

As the spins S a lie in the 2-3-plane within the frame S, 
we can write 

A^^A^OA,*,), (3.38) 

where a — 1,2. Notice that the angles 9 a are fixed, as 
the length of the edges of the triangle are given by Si :— 
\Si\, S 2 := |S 2 |, and S := |S|. We further introduce 
co-rotating frames C a for the spins S a . The remaining 
transition to the co-rotating frames is parametrized as 

KjJc. =A ij (0,f3 a , 7a ). (3.39) 
The decomposition Ay c = Aj S A s ^ leads to 

n) Ca = Ay s n SCa + n) s , (3.40) 

cf. (3.26). Then the kinematic terms can be written as 

E \SaiMk = S^) S + E S a^SC a ■ (3-41) 
a a 

These terms can be evaluated as in Sec. IIIB3. For the 
first term it immediately follows that 

fifaljg = Sa s + S cos 9 s <f> s ■ (3.42) 



7 



In the other terms we insert 

SI = S.A&A|| a , (3.43) 

*ssAc a =rts.c.-rts.s> (3-44) 
and finally obtain 

E 2 S ™&W a = ^ ^ + Sa S + ^(Sacia + S a a a ) , 
a a 

(3.45) 

where S 3 :— S cos 6$ and S 3 C := S a cos/3 a . The Poisson 
brackets read 

1 = {0 S , S 3 } = {a s , S} = { 7a , S 3 aC } = {a a , S a } , 

(3.46) 

all other zero. The dimension of the phase space is 12, 
as expected. 

Remember that the angles 9 a are not part of the phase 
space, but must be fixed from geometrical considerations 
in terms of the other variables. As the vectors Si, S2, 
and — S form a triangle and its lengths are part of the 
phase space, the angles 9 a can be obtained from the law 
of Cosines (see also Sec. IIIC2). 

C. Canonical angles: binary system 

We are now going to construct the phase space of a 
binary system with spinning components. We have just 
seen how the spins Si and S2 of the components combine 
to the total spin S := Si + S2. Now we also add the 
orbital angular momentum L to form the total angular 
momentum J of the binary system, 

J = L + S. (3.47) 

Notice that J is conserved to the post-Newtonian order 
considered here, which we use to reduce the number of 
variables. 

Notice that one could build the phase space simply by 
(3.36) and two copies of (3.32). However, as we want to 
make use of the conservation of J later, it is convenient 
to include J in the phase space. 

1. Complete phase space 
In Sec. IIIB5 we found that 

E \ S ^c a = + E^7a + S a a a ) . (3.48) 

a a 

We will now apply the same procedure to the terms 

L l n] L + S l n] s . (3.49) 
That is, we introduce a frame J by 

A# = A«(^, 0.,, aj), (3.50) 



such that J 1 — JAjj (here J := |J|) and the vectors L 
and S are lying in the 2-3-plane in the frame J. The 
parameterization of the frames L and S is now given 
relative to the frame J as 

A j j, = A«(0, 6 L , a L ) , A*/ s = A«(0, S , a s ) . (3.51) 

The result from Sec. Ill B 5 translates into 

VVL\ L + S l {l} s = J 3 <j)j + Jaj + La L + Sa s , (3.52) 

where J 3 :— J cos 9j. The angles 9l and 9s must be 
obtained from geometric considerations. 
In total we have 

PiT 1 + E \ S a-ij n W a = PrT + J 3 <j>J + 5?c7l + Sfc 72 
a 

+ Jctj + La L + Sa s + Sia 1 + S 2 a 2 , (3.53) 
and the Poisson brackets read 

1 = {r, Pr } = {<j)j, J 3 } = { 7 iA 3 c} = {12,S 3 C } 
= {aj, J} = {oil, L} = {a s , S} (3.54) 
= {ai,Si} = {a 2 , S 2 } , 

all other zero. 



2. Reduced phase space 

As in Sec. Ill A we are now going to utilize conserved 
quantities to reduce the number of phase space variables. 

We assume that the Hamiltonian is independent of the 
orientation of the co-rotating frames of the two com- 
ponents. That is, the Hamiltonian is independent of 
7a, S 3 C , and a a . Then the canonical conjugates S 3 C , 
7a, and S a are constant. The corresponding kinematic 
terms therefore turn into total time derivatives and can 
be dropped. 

Next we make use of the conservation of total angular 
momentum J. This allows one to align J with the 3-axis, 
i.e., 

J 3 = J & 9 j = . (3.55) 

In order to make contact with the notation in [55] , intro- 
duce the alternative notation 

aj + 0j = T , a L = Lp , a s = (j) S - ~ , (3.56) 

6l = 6 , 9 S = 6 - a ks , 9 X = s , 9 2 = s + a 12 - tt . 

(3.57) 

The kinematic terms simplify to 

p t r l + E l S «ijtt% a = VrT + JT + Llp + S<P S , (3.58) 

a 

and the Poisson brackets are 

1 = {r,Pr} = {T, J} = {<p, L} = {05, S} , (3.59) 



all other zero. There is still one conserved quantity left, 
namely J. This means that the Hamiltonian is indepen- 
dent of T, and one could even drop J and T from the 
phase space, too. However, T is needed for gravitational 
wave forms, so we will keep it for now. (Similarly, if one 
could observe the absolute orientation of the components 
somehow, then one would like to keep the variables a\ 
and Q!2, too.) 

Notice that beforehand we have formulated the phase 
space in a way that optimally allowed us to implement 
the conservation laws as constraints. This corresponds 
to the change from Lagrange I to Lagrange II in such 
a way that we have found generalized coordinates that 
eliminate constraints on the dynamical variables. 

Let us summarize the transition to the new variables, 
which must be inserted into the Hamiltonian: 

r' =rA ,:l (T,e,</?), (3.60) 

Pl = p r A a (T, 9, <p) + -A l2 (T, 8, tp) , (3.61) 
r 

S{ =S 1 A»(T, 6 Sl ct> s -K/2)V z (0,^,0), (3.62) 

S3 - S 2 A«(T, 6 s ,4>s- ^/2)A^' 3 (0, 9 2 , 0) . (3.63) 

Remember that we need to solve for 6, 0$, #i, and 2 , 
using geometric considerations (we know the length of all 
edges in the triangle where they appear, so one can apply 
the law of Cosines and Sines) , 



cos 6 = 
cos 9 s = 
cos#i = 

COS 02 = 

Si sin 6*i = 
L sin 9 



J 2 


+ L 2 - 


S 2 




2JL 




J 2 


+ S 2 - 


L 2 




2JS 




s 2 


+ s 2 - 


si 




2SSi 




s 2 


+ sf- 


s\ 



2SS 2 
-S 2 sin 2 , 
-S sin S . 



(3.64) 
(3.65) 
(3.66) 

(3.67) 

(3.68) 
(3.69) 



Notice that from now on we are not utilizing upper in- 
dices to denote vector components any more, upper in- 
dices are always exponents. Alternatively for s, a^ s and 
a\2 we have, see [55] and Fig. 1, 



«12 = cos 




(s 2 + s 2 -s 2 )M 

1 a n5 n5 ' 



4S 2 S 2 



/ Q2 Q2 _ Q2 

— 2S1S2 



■ -1 j L (J 2 + i 2 -s 2 ) 

= 7r - sm 1 — \ / 1 - — — r 



s 



4J 2 L 2 



(3.70) 
(3.71) 
(3.72) 



However, we will actually not utilize these inverse 
trigonometric functions, as they only give unique val- 
ues if additional assumptions on the spin orientations are 



made. Still Poisson brackets can be uniquely evaluated 
as follows. Partial derivatives of the constrained angles 
with respect to the canonical variables can be obtained 
by differentiating the law of cosines, e.g., 



sin9(J,L.S) 



dQ( J, L, S) 
dS 



S 

Jl 



(3.73) 



We actually substitute the Cosines of the constrained 
angles in terms of canonical variables using the law of 
Cosines. The Sines are first reduced to sinews and sinai2 
using 



(3.74) 
(3.75) 



sin 


sin s 


sin a ks 


S 


L 


J 


sin 0i 


sin 02 


sinai2 


s 2 


Si 


S 



Then it holds 

<9sinai 2 (S, Si, S 2 ) 
dS 

9sin«ks('/, L, S) 
dJ 

d sin a^ s {J, L, S) 
dL 

d sin «ks(</, L, S) 
dS 



S(S 2 + S 2 2 - S 2 ) 
2S 2 S 2 sinai 2 (S, S X ,S 2 ) 

J(L 2 + S 2 - J 2 ) 
2L 2 S 2 sina ks (J, J L,S) ' 

(J 2 - S 2 ) 2 - L 2 
4L 3 S 2 sin a ks (J,L,S) ' 

(J 2 - L 2 ) 2 - S 2 
4L 2 S 3 sin a ks ( J, L,S) ' 



(3.76) 
(3.77) 
(3.78) 
(3.79) 



These relations are used extensively throughout the 
present publication for the calculation of Poisson brack- 
ets. 




L 



FIG. 1. The geometry of the rotation angles. The dashed line 
represents the vector from the center of mass to the reduced 
mass. The vectors (i, j) span the orbital plane perpendicular 
to := j^. The vectors (i,js) span the plane perpendicular 
to k s := f . 
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IV. THE METHOD FOR ELIMINATING 
PERIODIC TERMS IN CANONICAL EOM 

In this section, we summarize some main ideas how to 
eliminate purely periodic perturbation terms in Hamilto- 
nian functions explained in [58, 59]. The linchpin there is 
to apply finite Lie transformations with a special choice 
of the generating function. 

Let the set X a represent the canonical coordinates and 
momenta, say x a — X a and p a = X n+a for a = l...n. 
Having chosen a generating function s, the transformed 
quantity X will follow as 

X = X + {X > s} + ±{{X,s},s} + ... 

= X + v s x + ^V 2 S X + ... 



+ 2^ki' D sX^.e^X, (4.1) 

k>l 

V°X := X , (4.2) 

VlX = {X,s}, (4.3) 

V 2 S X = {{X, S }, S } , (4.4) 



and so on. The Poisson brackets {X a , X n+a \ turn out to 
be invariant under this transformation, see, e.g., [58, 61]. 
It is therefore just a particular canonical transformation. 

To make this transformation valuable, we let the 
Hamiltonian function consist of an integrable part F* 
and a term i?W which is purely oscillatory (as the solu- 
tion to F* is inserted), 

HpQ = F*(X)+R^(X). (4.5) 

with O(iiW) = e and e < 1 as well. For the dy- 
namics of F* alone, an analytical solution to the sys- 
tem of canonical equations of motion may be known, 
X a . solution — : X a (t). Then, via a contact transforma- 
tion, we are able to shift the perturbation term to order 
C(e 2 ). This works well if we choose the generating func- 
tion to be 

s= I R\x=x(t>)dt' (4.6) 

J to 

and re-express the result purely in terms of canonical 
variables again. That is, the explicit time dependence is 
removed with the help of the solution X(t') = X after 
the integration was performed with the help of inverting 
the relation 

X = X(t + Cl ,C 2 ,...,C 2 n), (4.7) 

with the c's as a set of 2n integration constants, to finally 
get 

t + ci = ip x {X), (4.8) 
c j = <p j (X),(j = 2...2n). (4.9) 



From Eq. (4.6) one obtains 

{F*,s} + R = 0. (4.10) 

Note that (4.6) holds as we employ a canonical supple- 
mentary system where only the integrable part of the 
Hamiltonian fixes the evolution. After re-expressing each 
component of the generator s in terms of canonical coor- 
dinates only, it does not matter using which dynamics s 
has been computed, - it has the nice feature to definitely 
eliminate the leading-order (in e) terms of the perturbing 
Hamiltonian, see the computation below 2 . Performing 
the integral, the new Hamiltonian H written in terms of 
the transformed variables X will have the form (taking 
the theorem of interchange-of-variables in the Hamilto- 
nian as input) 

H(l) sp i H R(X) +Y -*-V k s (F* + R) 

k>l 

splitk=l + R ^ + p i F * 

V v ' 

=R+{F* ,s}=0 
k>2 ' fe>l 

k>i 

• > 

— : K, now of order C(e 2 ) 

because O(R) = e and O(s) = e. In other words: the 
new Hamiltonian will be just the old integrable part (with 
simply the new variables as arguments) and a higher- 
order series involving Lie derivatives of the oscillatory 
part. Let us call this new term K (X). This last part may 
contain terms which are integrable and can be absorbed, 
for a second transformation, into F*^ as 

F*^(X) := F*(X) + integrable part of K{X) . (4.12) 

This may be repeated until the remaining error is only 
of some required order. For the reader's convenience, 
an application of the method to the Duffing-oscillator to 
fourth order in the smallness parameter is provided in 
[59]. 

A remark on the form of our rest term: The non- 
integrable remaining term that appeared in Ref. [59], see 
their Eq. (5.33), has been modified in our computation. 



2 It is a vital detail that _R is assumed to be oscillatory only, see 
[58]. 
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The reason is the appearance of terms of higher orders 
in e - emanating from their definition of s - that we are 
unable to control. 

A further remark on our notation: The coordinates 
X a and the Hamiltonians F* will get an ordering num- 
ber - superscript (n) - after the n th transformation. The 
terms R will have a superscript (1) from the beginning 
on to emphasize that they already contribute to the first 
transformation and are at least of the lowest order of the 
smallness parameter. All (say n) performed Lie trans- 



formations have to be combined to obtain a connection 
of the final JfW to the initial X<-°\ This can be ob- 
tained by Eq. (5.14) of Ref. [59], where it is stated how 
to combine two Lie transformations with two generating 
functions to get only one, and which can be extended to 
multiple transformations. 

For the case of spinning compact binaries evolving in 
circular orbits, the problem of solving the equations of 
motion linear in spin will turn out to fulfill our require- 
ments to work with Lie transformations to get rid of com- 
plicated oscillatory terms. 



V. APPLICATION TO A SYSTEM WITH TWO SPINS 



A. Involved Hamiltonians 



The point-mass Hamiltonians to second post-Newtonian accuracy [62] will be given below. 



ttN Pr t L 1 

HpM - T + 2^ r 



H 



1PN 

PM 



2 fL 4 (3?7-l) £ 2 (t/ + 3) 1 ; 2 / (3t/ - 1)L 2 2t/ + 3\ 1 



8r 4 



2r 3 



2r 2 



Pr 



4r 2 



2r 



+ „(3t?-1)p: 



r f 



tt: 
xl 



2PN 
PM 



_ 4 |i 6 (5(r?-l)?7+l) L 4 (3t/ 2 + 20t/ - 5) | L 2 (8t/ + 5) , -3ry - 1 



16r 6 



8r 5 



2r 



4 



4 r 3 



1 



4 / -8?/ 2 - 20r/ + 5 3(5(r; - 1)77 + 1)L 2 



T6 V " V ' ~ rl ' ' ' rr V 8r 16r 2 

/ (— 47/ 2 - 20?/ + 5) L 2 3(5(7/ - 1)7/ + 1)L 4 11t/ + 5 X 



+ Prl 4r 3 1 16r 4 

The terms linear in spin through next-to-leading order read [9] 



2r 2 



ttSO 



ttSO 
n NLO 



4 r 3 



{(277 + 3^1-477 + 3) (L ■ Si) + (277 — 3\/l — 477 + 3) (L-S 2 )} 



16r 4 



(L-Sj 



r ((6t/ 2 + 26t/ - 5) p 2 + 677(277 + l)p 2 r ) 



+ y/l - 4t/ (r (67/p 2 + (16t/ - 5)p 2 ) - 8(2t/ + 5)) - 8(4t/ + 5) 



-(L-S 2 ) 



r ((6t/ 2 + 267/ - 5) p 2 + 67/(277 + l)p 2 ) 



+a/1 - 4t/ (r ((5 - 16t/) P 2 - 6t/p 2 ) + 8(27/ + 5)) - 8(4t/ + 5) 



(5.1) 
(5.2) 



(5.3) 



(5.4) 



(5.5) 



where p 2 = (n i2 ■ p) 2 + (n J2 xp) 2 =pj+^, making use 
of p r := (ni2-p), r = ni 2 r, and L := r x p. None of these 
terms depend on the orbital tp and the magnitude L is 
unaltered. The spin-orbit Hamiltonians can be expressed 
in terms of 

S := Si + S 2 and A := Si - S 2 , (5.6) 
where one can represent the scalar products (L • S) and 



I 

(L • A) in terms of the new canonical variables. 
This gives 

(L • S) = \ (J 2 — L 2 — S 2 ) , (5.7) 
(L • A) = ^{(Si - S 2 )(S 1 + S 2 ) (-L 2 + J 2 - S 2 ) 

-wa.(<l>s)A{J,L,S)}, (5.8) 
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having introduced the shorthand 



A{ J, L, S) — ALSSiS2sina^ s (J, L, S) sm.ai2(S, Si, S2) , 
=:A. (5.9) 



For further computations we also introduce the abbrevi- 
ations 



Gi{J, L, Si, S2) :— {Si — S2) (Si + S2) 

x ( J — L)(L + J) . (5.10) 



and 

Aj = -4L 2 S 2 sin 2 a ks ( J, L, S) , (5.11) 
A s = -AS 2 S 2 sin 2 a 12 (S, Si, S 2 ) . (5.12) 

It holds A 2 = A j A s and 

Aj = {J-L-S){J + L-S) 

x(J-L + S)(J + L + S), (5.13) 
A S = (S-S 1 -S 2 )(S + S 1 -S 2 ) 

x(S -S 1 + S 2 ){S + Si + S 2 ). (5.14) 

The following subsection will deal with circular orbits and 
the question how one can evaluate the Poisson brackets 
imposing this restriction. 



B. Evaluating the Poisson Brackets in the Circular Case 

For evaluating the Poisson brackets, we find it important to mention that for some variable 
X G M. = {ip, L, Y, J, 4>s, S, a a , S a } (M. is the phase space elements without p r and r) and an arbitrary 
function Z of Ai it holds 



{X,Z} 



r— f(L,...) 



p.r=0 

r—r(L....) 



(5.15) 



which means that we can evaluate {~,^} for two quantities on a general orbit and later impose the condition of 
circularity, replacing r by its solution f(L, ...) to p r = 0, or we can do it the other way round. The case of circularity 
is no primary constraint and works without artificial/external forces. The solution to p r = can be given as follows: 

L 2 - 4c" 2 

21 ( J 2 — S 2 ) 21 / 9(Asin(0 s )- (J 2 -S 2 ) (S 2 - S 2 )) 9 (Sj - S 2 ) \ 
16 + 6 \ 4L 2 S 2 4S 2 J 

43e 2 253 \ <c _J 1605 t 1605 ( J 2 - S 2 ) 





8L 2 32L 2 
J 2 -S 2 )(S 2 -S 2 )) (S 2 -S 2 ) 



L 4 S 2 L 2 S 2 

(Asin(^)-(j 2 -S 2 )(S 2 -S 2 )) t (S 2 -S 2 )' 





L 2 S 2 



(5.16) 



As we decide for the second choice in Eq. (5.15), this is inserted into the total Hamiltonian H, which will also generate 
terms of order 0(c~ 3 ) and 0(c -5 ). The Lie transformation terms and the terms dictating the angular velocity ft s 
(defined by 4>s = 0, s t + <j)so) will origin only in those odd powers of c. We next set 

V=--\-e 2 (5-17) 

where e <C 1 in order to emphasize that we work in the region of almost equal masses. We are able to write down the 
circular-orbit Hamiltonians that contribute to the spin-orbit part (=those having odd powers on c) in terms of e, 

pSO.circ -Jl (T 2 T 2 2 ^ 3 £ (S?-Sf) e 2 7\ 3eAM4>8) \ 

= c )T^ J - L - S H AS 2 4 16 " 4L*S 2 > ' (5 - 18) 
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w so.circ_ c- 5 j { T 2 T 2 C 2\ f 99e (Sf - SI) 3e 4 93e 2 975 \ 99eAsin(0 s ) | 

Hnlo -^^ J - L - S >y ^52 + T6^^2" + 256^1 16S~ 2 /' (5 " 19) 

and for the sake of completeness, those coming from the point-mass parts only, 



H PM,c ir c = __ j (5 20) 



1 



jPM, circ — 2 



w PM,circ_ _ 4 -16e 4 -104e 2 -1269 

H 2PN -c , (5.22) 

The expressions (5.18)-(5.22) are not a Taylor series in e where higher orders have been neglected. In fact, higher 
orders do not exist. 

C. Initial Decomposition 

We shall sketch the Lie transformation procedure representatively for the first order and truncate our procedure 
to the order C(e 3 ), because terms of order C(e 4 ) are in direct competition with the second-order Lie transformation. 
The first "perturbing" Hamiltonian with purely oscillatory character may be all that comprises terms of HgQ C with 
sin(0s): 

77!*(0l TT PN,circ . TT lPN,circ . TT 2PN,circ 
r = % + H N + H N 

?2 



L 6 



c 



( M %# ll£ + ^(7"4e 2 ))}(J 2 -i 2 -5 2 ) 

248e 2 + 325) ) 1 ( J 2 - L 2 - S 2 ) , (5.23) 



f99(S 2 -S 2 )e , 3 , ~ 



L 8 \\ 16S 2 256 
3/1 

R {1) = ~YeiJ^ sin(0s) i 41 ^ 6 + 33c ~ 5e ) > ( 5 - 24 ) 

because, as we compute the Poisson brackets of S and 05, we see that 

j 5;F *(o)j = 0; (5.25) 

{^,F* (0) } = ^ 0) , (5.26) 

O ( 0) := g 1 1 s (4e2 _ 7) _ ^(j.L.gx.g,) 

■ ',;{ ^^;^> ,;,( 2„ 2 .325)} 

= const. (5.27) 

in (applying our sense of rotation) accordance with the equal-mass case (e — > 0) evolution equation for (f>s, Eq. (4.6d) 
in Ref. [55], where geometrical considerations lead to 

L/sin(a ks - 6) = 5/ sin 6, (5.28) 

see Fig. 2. The first generating function is defined to be the time integral of to some time t, so we set tfis — 
Q^t + cfrsQ with some irrelevant 4>so an( i S = Sq; then we perform the time integral, and re-express Q^pt + <f>so ~* 4>s- 
Next, we compute the new Hamilton function according to Eq. (4.11) and get F*W(JfW) = F*W (X^) + (X (1) ) 
with 



K <2 >(X (1| ) = i {«,*<"} 



2 
3! 
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FIG. 2. The plane spanned by L and J. Adding auxiliary lines, we immediately see that (5.28) holds. 



+ ... (5.29) 

Notice that the R above possesses integrable contents. These contents are all that remains when one removes the 
trigonometric functions of 4>S m 

= a(X^) /3 m (X^)sm(m(f>s) +jn(X {1) ) cos(ncfe) , {m,»i}eN. (5.30) 

mtcgrablc purely periodic 

As we absorb these terms into F*V-\ the new Hamiltonian after the first transformation will look as what follows 
after a short break where we compare our calculation to the aligned-spin case for convenience of the reader. 



D. Comparison to the aligned-spin case 



To strengthen the faith in our system of canonical coordinates and conjugated momenta, and to using the Poisson 
brackets after imposing circular orbits, we may do a small check. In the case of circular orbits at "almost alignment" , 
we see that the spherical coordinates (r,<fi,8) satisfy r — const, 8 < 1, 9 w |, and <fi = T + tp + 0(Q 2 ), see Eqs. (4.27) 
of [52]. As well, we can impose S ~ J — L and 9< 1, which can be justified by conservation of J and what triangular 
relations for the polygon (J, L, Si, S2) constrain, for further calculations. 

In [63] , it has been assumed that 777,1 > 7712 , which is in contradiction to our assumption. This changes two signs in 
the function 77 = 77(7711,7712) compared to our computation. We also have to re-scale the spins according to the rules 
given therein. Note that the conservation of the dimensionless J is only possible if the individual spins are scaled in 
the same way L is. 

Having done this, we compare our result (TSS) for the spherical phase velocity to the orbital phase velocity taken 
from [63] (THS), defining <I> to be the elapsed total phase through one radial period V r , and see - through 1PN order 
linear in spin - 

* - 1 ■ W ^ + 5( gl+M _ 10e(gi _ ga)+ 37L-35( gl + ga )\ (531) 



^THS.circ LS C 2 V 2 



ll'C— sc 



1 ; 1 / e 2 (L + 5(S*i + S 2 )) in 37L - 35(5, + S-. 



"TSS,eire ~ (J + ^ = + ( -^^^^A IQ^ - <? 2 ) + — ' ^ ) , (5.32) 



which shows agreement. 
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E. After the First Transformation 



To omit the superscript "(1)", the arguments X a are we replace each term by its transformed version, X a 
on both sides of the equation: 

p*(i) = p*(0) 



+ c- 3 <Aj[ A s 




243e 3 G 4 (J,L,Si,S b ) 45e 2 \ 



54 



98L e S w 



56L e S 6 ) e 49L 6 



c2 
°b 



1 



s s ' s s s e 



G 4 {J, L, Si, Sb 



i _s f_sl 
s 2 S 4 s 4 

1 J 2 



L G S 2 L 6 S 4 L 4 S 4 




1 



1 



Gi{J, L, Si, St,) 



L 6 S 8 L 6 S 6 L 4 S 8 

99873e 3 G 4 {J,L,Si,S b ) 39285e 2 \ 



5488 J L 8 S' 10 



6272L 8 S' 6 J 



. 11097 
137217 



1 



C2 c2 
1 r_ b _ 

s* 8 s 8 s e 



A, I e 3 



11097 
1372 



J 



1 



Gi{J,L,Si,S b )+( 2 
1 



7857 


' 1 


si 


c2 
^b 


1568L 8 


s 2 


s 4 


S 4 



L 8 S 8 L 8 S 6 L 6 S 8 



Gi{J, L, Si, Sb) + € 



,7857 
1568 



1 



J 



1 



L 8 S 2 L S S 4 L 6 S 4 



= (e 2 ) 



x„ 



iX 



(1) 

(5.33) 
(5.34) 



After the Second Transformation 



As one performed a second transformation, the perturbing function R would be (after the split, naturally) shifted 
to fourth order in e, and the integrable part also absorbed terms of order e 4 , 

p*(2) ( X (2)\ = F *(l) f x (2)^ 



R^ fx( 



(2) 



G(6 4 ) 



(5.35) 
(5.36) 



We will not perform this transformation and stop the calculation here. All emanating residues are of the form 
~ sin(m^>s) and ~ cos(m(j)s) with m and n as positive integers and having complicated functions of spin amplitudes 
as total prefactors, which one can easily verify for all orders. The residue after the first transformation is of order 
C(e 2 ) - when talking about (5.33), we speak of an integrable system of first order in the perturbation parameter e. 
Its solution (labeled with a "bar") reads 



{4>f,F*W} 



= {T( 1 ),F*( 1 )} , 



(5.37) 



(l> 8 (t) = n s t + <M* = o), r> s 
<p(t) = n L t + ip(t = o) , n L 
f(t) = n r t + r(t = o) , n T 
S{t) = S(t = 0). 

Because the generator only affects the remaining term, the integrable part is unaffected, and the circular point mass 
Hamiltonians always keep their form as they belong to F* (up to the fact that the variables get new names). The 
reader should note that there are "Newtonian" terms in the generating functions, clearly speaking: terms of order 
O(o ), which arise because of the fact that the precession velocity is of the order O(c~ 2 )O(e ), and when integrating 
R, these velocities become some part of the denominator. 



VI. COMBINING FURTHER CANONICAL TRANSFORMATIONS 



From Section 11.2.3 of Ref. [59] we know that combining two Lie transformations with generators and s^ 2 ' will 
be expressible as performing a single Lie transformation with the generator 



,(3,1) = ( S (D + s (2)) + 1 { S (D )S ( 2 )} + 1 | S (D _ s (2) j { S (D, S ( 2 )}} 



(6.1) 
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The above transformation connects the variables X^ 2 ) (those after the second transformation) to the initial ones X 
via 

x = x< 2 > + i {x^,^ 1 '} + i {{x( 2 > s w>} lS w)} . . . (6.2) 

The transformation may be inverted and the resulting s^ 2,1 -* be expressed entirely in terms of the initial X. We 
will skip this formula because of reasons of comprehensibility and state that, if the reader is interested in further 
transformations due to the reduction of oscillatory remainder functions, the total generating function of n successive 
Lie transformations can be obtained from the recursion scheme 

a™ = s« , (6.3) 
S ( 2 'D = ( S W + a W) + l - { S «, S ( 2 )} + 1 { S W _ s (2), { S W, S ( 2 )}} + ... , (6.4) 

s (»,n-l) = (s (n-l) + s (n)) + * | a (n-l) > fl (n) j + _L j>~D _ s («) ) j>"l) , 1 1 + ... . (6 . 5) 

Structurally, as we would take the full spin-orbit Hamiltonian as input to our scheme without truncating after third 
order of the smallness parameter e, the residues emanating after the n th iteration are going to appear at the following 
orders of e: 



step # 


c 1 


c 2 




c 4 


e 6 


e 6 


c 7 


e 8 ... 





* 




* 












1 




* 


* 


* 


* 


* 


* 


* 


2 








* 


* 


* 


* 


* 


3 
















* 



TABLE I. Position of terms contributing to the oscillatory residuum, starting from the untruncated spin-orbit Hamiltonian, 
Eqs. (5.18) and (5.19). Step means the initial form. From step 1 onwards, there are infinitely many terms at higher orders 
of e. 



VII. SOME REMARKS ABOUT HIGHER ORDERS 



In [58] it was stated that the Lie series converges if there exists a finite number B such that the generating function 
satisfies 



\s\ < B . 



gki+k 2 ...m 1 ...+r, 



^■'V 'V -'V 



< B 



(7.1) 



For simplicity, we use the perturbation function at linear order in e and first PN order, which is of structure 



3e (51c 



112L 2 ) Acos((j)s) 



392L 2 S 3 



0{e 3 



(7.2) 



However, the quantity A generates more and more terms at each evaluation of the Poisson bracket which grow in their 
magnitude. One can see this fact as well in Eqs. (3.76) - (3.79), where the sines can have values arbitrarily close to 
zero. This circumstance violates Eq. (7.1) and thus generates - in general - an asymptotic Lie series. It depends on 
the system (i.e.: the initial configuration) and on the mass parameters how many terms of this asymptotic series can 
be taken to properly describe the dynamics of the spin-orbit problem 3 . 



3 A first numerical insight showed that the range of S(t) due to the tQ low _ range oscil i ations around the initial spin length S(t = 0), 

initial Hamiltonian depends on the initial value <j> s {t = 0). Some which may algo aflfect the speed of divergence . 

values lead to the full range J—L < S < 51+52 while others lead 
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Let us give a numerical representation of the results and their speed of divergence through a sequence of 3 modified 
Lie transformations for an example set 

X num = Is = 1, Si = 1, S 2 = i L = 10, J = 111 . (7.3) 

Then, after 3 iterations, we obtain a value of 
p*(3) _ _q.0050 

num w.ww^w 

+ cr 2 l(T 4 {-1.15 + 0.13e 2 } 

+ c~ 3 {0.0000086 + 0.0000113e + 0.0001276e 2 - 0.0035799e 3 + 0.0377074e 4 + 3.2521845e 5 } 
+ C - 4 10~ 6 {-4.957 - 0.406e 2 - 0.063e 4 } 

+ c~ 5 {0.00000076 + 0.00000093e + 0.00000975e 2 - 0.00026274e 3 + 0.00259869e 4 + 0.20891229e 5 } 

+ 0(e 6 ) . (7.4) 



VIII. CONCLUSIONS AND OUTLOOK 



In this article, we have found a reduced variable space 
for the treatment of the binary spin-orbit interactions. 
This space consists of pairs of compact angle variables 
and their conjugate momenta, which are amplitudes of 
angular momenta defining the plane in which the specific 
angle is evolving. 

Taking all spin-orbit interactions through NLO and 2PN 
point mass contributions, we could solve the resulting 
equations of motion fir the circular-orbit case taking as 
help the Lie transformation method to successively get 
rid of (until now) untreated oscillatory terms. It showed 
off that, by construction, the structure of the point mass 
Hamiltonians is not affected by the Lie transformation 
algorithm as there are no (S Q • 1112) couplings in the spin- 
orbit terms. 

Outlook: A numerical comparison of the transformed 
Hamiltonian to the original one may give some insight 
about the correctness with respect to the order of the 
perturbation parameter e, regarding a number of con- 
figurations. We also concern resummation techniques in 
the future that keep the structure of singular points in 
the perturbing Hamiltonians but remarkably reduce the 
number of the involved terms. 

It should be straightforward to extend the results of 
the present paper to higher post-Newtonian orders, i.e. 



to include the NNLO spin-orbit Hamiltonian. It would 
be valuable to consider other methods and perturbation 
parameters e to analytically obtain solutions to the equa- 
tions of motion. As well, as our phase space generates 
an asymptotic series for the spin-orbit problem, a new 
proposal for canonical variables (where no sines of non- 
canonical orientation angles, but polynomials appear in 
the generator) may mean a future improvement. 
This can also facilitate an extension to higher orders in 
spin possible, e.g., to spin(l)-spin(2) interactions, which 
actually seem to be impossible to be handled with the 
current approach. Furthermore an extension to eccen- 
tric orbit by an expansion around circular-orbit case is 
envisaged. 
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